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Abstract 

We calculate the non-universal correction to the Zbb vertex in a simple ex- 
tension of the Standard Model, where a charge +2/3 isosinglet quark is added 
to the standard spectrum. Comparison is made with other solutions to Rb (and 
R c ) that demand particles lighter than My/- 



' Permanent address. 



Current experimental data on Rb (= T(Z — > bb)/T(Z — > hadrons)) and i? c deviate 
significantly from Standard Model (SM) predictions [|TJ. It has been suggested that both 
discrepancies can be accounted for if one makes a minimal extension of the SM by adding 
an extra up-type isosinglet quark Q §. In this scenario, R c is reduced because of c-Q 
mixing, while Rb is enhanced by a more elaborate mechanism. The top quark is light, and 
remains hidden below M\y (i.e. m t < Mw) by some fast decay mechanism, while the 180 
GeV quark observed at the Tevatron is identified as the dominantly singlet quark. As a 
result, the effective top quark mass, defined as the equivalent m t that appears in the Zbb 
vertex correction within the Standard Model, is reduced and thus Rb is increased. The main 
purpose of this paper is to perform a detailed calculation of the non-universal correction 
to the Zbb vertex when Q is present. The novel feature is the presence of tree level flavor 
changing neutral current (FCNC) Z-t-Q couplings in the loop. 

Adding an up-type isosinglet quark gives rise to new gauge invariant mass terms of the 
type M QIQr an d M'- Q L u® Rl where denotes the gauge eigenstates of standard up-type 
quarks. Moreover, there are also mass terms of the form m! i u° iL Q R which arise from additional 
Yukawa couplings. Since we do not expect the first generation to play an important role 
in the present context, we will ignore them throughout the paper. We therefore have the 
following mass matrix for the up- type quarks ||, 

which can be diagonalized by a biunitary transformation as usual, 

M = diag(m c , m t , m Q ) = ^M T U . (2) 

In terms of the mass eigenstates u and d, the charged current takes the form [|J, 

u a L (V aj ) (3) 

where a ranges over c, t, Q, and j = s, b. The quark mixing matrix V is now 3x2 and 
hence non-unitary. Since \V c b\ 2 — (0.04) 2 <^ 1, we shall set V c b = throughout the paper, 
which leads to considerable simplification of our results. 
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The neutral current for ti-type quarks now becomes ||], 

u L (U t%-I Sw) + u R (-1 s^) 7 m m r , (4) 



where we have suppressed flavor indices, U = VV\ and sw = sm9w is the effective elec- 
troweak mixing angle at the Z-pole. Only the isospin part is modified since only Q° L is 
nonstandard. Because V is non-unitary, U has off-diagonal terms. We also note that in 
general, \U CC \ < 1, and can be used to account for the apparent suppression of R c 0. 

To calculate the Zbb vertex correction, we must evaluate the ten Feynman diagrams listed 
in Fig. 1. We make use of the REDUCE code set up by Hou and Stuart ||, which is in turn 
based on the package LERG-I || . The results are expressed in terms of two scalar integrals, 

d n k 1 



B (p 2 ;mj,ml) = J 

ri I 2 2 2 2 2 2\ 



m 2 (k 2 + m 2 ) [{k + p) 2 + m 2 ,] ' 
d n k 1 



in 2 (k 2 + m\ ) [(k + p\) 2 + ml] [(k + p\ + p 2 ) 2 + ml] 

The original program was set up for a general study of FCNC decays of the possible fourth 
generation b' quark. It was used by Lynn and Stuart |7| to calculate the Zbb vertex correction 
in the context of SM. They find that the genuine FCNC vertex diagrams are well-behaved 
in the limit that the masses of the two external quarks become equal, and the resulting ex- 
pression can be used without further ado. On the other hand, diagrams containing fermion 
self-energies individually would have spurious divergences in this limit, since one has a prop- 
agator factor of l/{m 2 — m 2 ,). However, the sum of the four diagrams has a well defined 
limit, and its value can be determined using the L'Hospital rule. 

In the on-shell scheme, the UV divergence is taken care of by the counterterm [0], 

(-2^ 3 cos e w 4/16n 2 )B (0; M 2 Wl M 2 W ). (5) 

One could equivalently make a subtraction at q 2 = M|. We shall not go into the details 
of the renormalization program, since we are concerned with only the internal quark mass 
dependent contribution. However, it is easy to understand the origin of the counterterm of 
eq. (|5p. In ref. ||, it was checked in great detail using elementary Ward identities that all 
the divergent pieces cancel except for a left-over piece coming from diagram (c) , which takes 



the form (2ig 3 cos8wti/lQ7r 2 )B (q 2 ; M^, M^). In the FCNC case, the term is removed by 
the GIM mechanism, but in the flavor diagonal case such as Zbb vertex, subtraction and 
renormalization is needed. Thus the counterterm in eq. (g) is the most natural one. 

For diagrams (c)-(j) there is only one quark in the loop, either t or Q, and the flavor 
factors are \V t b\ 2 and |Vq&| 2 , respectively. One therefore finds the contribution 

\V tb \ 2 F L ^ n \m t ) + \V Qb \ 2 F L W(m Q ), n = c, d, . . . , j, (6) 

where F^ n '{rn) is the SM result for diagram n, in the notation of ref. f|. On the other hand, 
as can be seen from eq. the two internal quarks in diagrams (a) and (b) can be different, 
which is an interesting feature of introducing isosinglet quarks. For identical internal quarks, 
the flavor factors are VXU tt V tb and V^UqqVqi,, while if the two internal quarks are different, 
the flavor factors are V t lU t QV Qb and VQ b UQ t V tb . 

It is useful to establish some relations between these flavor factors. Using the fact that 
U = VV^ and V^V = I, one finds that UV = V. Since we take V cb = 0, this translates to 

UttV a = V tb -U tQ V Qb , (7) 
UqqVqi, = V Qb — UQ t V tb , (8) 

in component form. Multiplying by V t l and VA b , respectively, from the hermitian nature of 
U, one finds that VQ b UQ t V tb and V^ b U t QVQ b are not only real, but equal to each other, hence, 

v* h u tt v th = \v tb \ 2 -v; b u tQ v Qb , (9) 
v$ b u QQ v Qb = \v Qb \ 2 - v; b u tQ v Qb . (io) 

We therefore have the following contribution from diagrams (a) and (b), 

\V tb \ 2 F L ^\m t ) + \V Qb \ 2 F L W(m Q ) 
- V; b U tQ V Qb {Fi n \m t ,m t )-2Fi n \m t ,m Q )+Fi n) (m Q ,m Q )} 
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According to eq. (^), only the t% part should be used for the calculation of the last term. 
Putting everything together, we finally obtain 

F L = \V tb \ 2 F L SM {m t ) + \V Qb \ 2 F L SM {m Q ) - V*U tQ V Qb £ AF| n) K,m Q ), (12) 



where 



AF^\m t ,m Q ) = {^ n) (m t , mt ) - 2FP{m u m Q ) + F^ n \m Q , 



Since V is part of a 3 x 3 unitary matrix, it can be parameterized as 

/ C 2 \ 



V 



— S2S3 Cj(6 l 



(13) 



(14) 



where Si = sin^j, Cj = cos6*j, and S2, S3 are the c-Q, t-Q mixing angles, respectively. 



Consequently, we have \V t 



tb\ 



Cl \Vc 



Qb\ 



I 2 = Sf , and V* b U tQ V Qb = CfCfSf. Note that the 
phase 5 is removable in any case because of setting V c b = 0. 

Inspecting eq. ( |1"2"D , we see that the first two terms are identical to the case of adding 



a fourth generation, with Q = t' and S3 = \Vfb\- F L is nothing but the full SM result of 
ref . Jf| , expressed in terms of three "universal functions" , 



SM 
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32n 2 c w 



{(t u 3 -Q u s 2 w )p-2tjc 2 v A-2t d 3 E} 



(15) 



The last term appears only for Zbb vertex but not in 766. We note that the expressions given 
in ref. [0] contain three typos, two of which can be identified simply by dimensionality. For 
m not much smaller than M\y, the leading effect is Fi SM {m) oc m 2 /My V . 

The AFl term is specific to adding singlet quarks, and arises from diagrams (a) and (b) 
only. We find that AF^ contributes only to the S term. This is hardly surprising since 
adding the singlet quark Q should affect only the Zbb vertex. The explicit result is 
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2 -a 2 
w Q 



m 2 Q ) C (q 2 ,0,0,m 2 Q ,m 2 Q ,M^) 



+- (M 2 ? -q 2 - ml) [M 2 W -q 2 - m 2 ) C (q 2 , 0, 0, m 2 Q , m 2 , M 2 W ) 

~ (M 2 , - q 2 - m 2 ) 2 C (q 2 , 0, 0, m 2 , m 2 , 
1 



+— {2M 2 W - 3q 2 - 2m 2 ,) [B (q 2 , m 2 Q , m 2 Q ) - B (q 2 , m 2 Q , m 2 ) 



2q 
1 



+— (2M 2 V - 3q 2 - 2m 2 ) [B (q 2 , m 2 , m 2 ) - B (q 2 , m 2 Q , m 2 



AE® 



rn, 



2M, 



2 C tf,0,0,ml,m 2 Q ,M$y)- 



m 2 Q m 2 t 



w 

4 

-^-C^MmlmlM^). 



C ( 



The leading effect comes from AH^, which takes on a rather simple form. This can be easily 
understood as follows. From diagram (b) with only part of the Z-a-f3 vertex, 

d 4 k itlR (/+ V + m a ) {j+ ft + m p ) 7 l 



(2tt) 4 (fc2 + M&) [(k + p'Y + ml] [(k + pf + m 2 
m\m% r d A k i« 7/l7£ 



MJv J (27T) 4 ( fc 2 + M 2^)[( fc+p /)2 +m 2] ( fc+p )2 +m 2 



2 2 

oc !^C (0,g 2 ,0;M^,m2,m5). (16) 



The sign of the AFl term can be understood as follows. As noted earlier, the first two terms 
of eq. flT2]) are identical to the case of adding a fourth generation, hence the AFl term 



represents the difference between singlet and fourth generation results. As Fl is mainly a 
measure of nondecoupling, one expects the AFl term to soften the effect, since in the singlet 
case the gauge invariant masses M' and M of eq. (1) are of decoupling nature, in contrast to 
the Yukawa masses m and m' . Note that AFl has a built-in cancelation mechanism, hence it 
is unlikely to be dominant. In Fig. 2 we compare F^ M (mt) and AF^(mt, toq) vs. mt for rriQ 
fixed at 180 GeV. Since 5| should not be too large, the AFl term is indeed subdominant. 
The non-universal correction to the Zbb width is [[| 

r(z-&&) = r°[i + ,f 6 ], (17) 

where T° includes all effects other than the one specific to the Zbb vertex. As a first approx- 
imation, taking leading effect only, our result is equivalent to an effective top mass 
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e ^ 2 ~ CM + 5 3 2 m 2 Q = to 2 + S 3 2 K - m 2 ), (18) 



in the SM Zbb loop, i.e. 5f M (m^ ), which can be fitted by 

e//\2 



<5f M (mp) = 0.01 -0.49^/- - 0.45 ] , (19) 



I TO, 

in the mass range of interest. The constant term differs from that of ref. || since it is 
scheme dependent. In Fig. 3 we compare our full result and the approximate one of eq. 
(Tl9|). The full result is smaller in absolute value by about 6 to 14 percent for given Sf value, 
in agreement with qualitative arguments given earlier. 



Since eq. (|T8| ) is different from what was used in ref. [Q, it is worthwhile to repeat the 
numerical estimates. S 2 is not affected and has the value 0.03 from taking the R c value of 
0.161. To fit R b 0.2219 with m t , m Q = 70, 180 GeV, but using eqs. (U) and flU), we 
find that m| ~ 95 GeV and S 2 ~ 0.15. The numerical value of S3, however, depends quite 
sensitively on mf , and hence on the small sub leading corrections that we have computed. 
Using the full result, we find from inspection of Fig. 3 that the actual S 2 is around 0.22. 
These numbers are not drastically different from ref. 0. 

It is argued in ref. that in order to hide the t quark and at same time keep Ri 
untouched, one has to introduce an additional Higgs doublet. As a result, there are additional 
corrections to the Zbb vertex from diagrams (b), (d), (i) and (j) due to the physical charged 
Higgs boson. Since the mass of the physical charged Higgs is taken to be much greater || 
than M\y, we expect its correction to the Zbb vertex to be suppressed. Indeed, by taking 
the physical charged Higgs mass to be 250 GeV, we have checked numerically that it only 
adds about 6-11 percent to the previous result of 5 b , since cot (3 is constrained by 5-mixing. 

We turn to some comments and discussion. The scenario of ref. [|| relates the Rb and 
R c problem to the existence of heavy singlet quark Q and light top quark. It is similar to 
the partial SUSY solution to R b (but not R c ) in predicting a host of scalars and fermions 
below My/. It is, however, in spirit closer to fourth generation models P,ID|. Note that for 
the fourth generation case, assuming that \V tb \ 2 + \V t , b \ 2 = 1, one has 

F L = F L SM {m t ) + \V t , b \ 2 [F L SM (m t i] - F L SM {m t )\. (20) 

Since the absolute value of FL SM (m) increases as m 2 , taking m t ~ 180 GeV and m t i > m t 
would aggravate the R b problem. To alleviate or resolve the R b problem, one must have 



either m t = 180 GeV and m t > < m t 0, or m t i = 180 GeV and m t < m t > ||10|| , in the 
convention that t is the isospin partner of the b quark. The models of refs. and |T(| also 
contain light supersymmetric particles, such as stop t and chargino/neutralino x ± , x°, which 
are needed either to hide the lighter quark, or explain the R b problem. 

Thus, in all these cases, light particles are predicted to exist below Mw, and should be 
readily discovered as LEP 1.6 turns on and accumulates sufficient integrated luminosity. As 
we eagerly await imminent discovery, we should keep in mind that all these scenarios could 
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be ruled out before the end of 1996. The implication would be that the Rb problem cannot 
be explained in the context of these models. It should be stressed, however, that the singlet 
model stands out in its ability to suppress R c with relative ease, without necessarily touching 
Rb- Experimentally, R c is harder to measure, hence it would take a long time before one 
can be confident that the experimental value is in full agreement with SM. Assuming that 
no light particles are found at LEP 1.6 and beyond, we could take mq > m t = 180 GeV 
and fit any small deviation in R c . If we take Sf ~ S| < 0.03, and/or if uiq ~ m t , we see 
from eq. (|T2|) that Rb is not much affected. It is amusing that if mq is not very different 



from m t , it could partially explain the larger ll a t f observed by CDF ||11|| . We urge the LEP 
experiments to continue the refinement of R c (and of course Rb) measurement. 

In this letter, we present a calculation of the non-universal correction to the Zbb vertex 
when a charge +2/3 isosinglet quark is added to the Standard Model. Since the GIM 
mechanism is violated, it is possible to have flavor changing neutral current Z-t-Q couplings 
in the loop. The result is close to but slightly weaker than that of adding a fourth generation, 
and can be approximated as an effective top mass. If we identify the dominantly singlet quark 
to be the one observed at the Tevatron and assume the top quark to be lighter than Mw, we 
can fit the the current experimental data of Rb (and R c ) by choosing appropriate S% (and 
Sf) values. If LEP 1.6 does not find the light top quark, the model cannot explain Rb but 
could still account for suppression of R c , even if mq > 180 GeV. 
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Fig. 1 The ten Feynman diagrams that contribute to the Zbb vertex. 
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12 3 4 

m t 2 /M| 

Fig. 2 Comparison of F^ M (m t ) (solid) and AF^\m t ,niQ), J2n=a ^Fj*\m t , mg), (dots, 
dash) vs. m t for uiq = 180 GeV. 
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Fig. 3 Comparison between the full result (solid) and the approximate one (dash). 5| is 
fixed at 0.03, while S| varies between 0.1 and 0.3, with m t and mq fixed at 70, 180 GeV, 
respectively. 
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